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ABSTRACT 

We calculate tlie eccentricity excitation of asteroids produced by the sweeping 
secular resonance during the epoch of planetesimal-driven giant planet migra- 
tion in the early history of the solar system. We derive analytical expressions 
for the magnitude of the eccentricity change and its dependence on the sweep 
rate and on planetary parameters; the sweeping leads to either an increase or 
a decrease of eccentricity depending on an asteroid's initial orbit. Based on the 
slowest rate of z/g sweeping that allows a remnant asteroid belt to survive, we de- 
rive a lower limit on Saturn's migration speed of ~ 0.15 AU My~^ during the era 
that the resonance swept through the inner asteroid belt (semimajor axis range 
2.1-2.8 AU). This rate limit is for Saturn's current eccentricity, and scales with 
the square of Saturn's eccentricity; the limit on Saturn's migration rate could 
be lower if Saturn's eccentricity were lower during its migration. Applied to an 
ensemble of fictitious asteroids, our calculations show that a prior single-peaked 
distribution of asteroid eccentricities would be transformed into a double-peaked 
distribution due to the sweeping of the vq. Examination of the orbital data of 
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main belt asteroids reveals that the proper eccentricities of the known bright 
[H < 10.8) asteroids may be consistent with a double-peaked distribution. If 
so, our theoretical analysis then yields two possible solutions for the migration 
rate of Saturn and for the dynamical states of the pre- migration asteroid belt: 
a dynamically cold state (single-peaked eccentricity distribution with mean of 
~ 0.05) linked with Saturn's migration speed ~ 4 AU My~^, or a dynamically 
hot state (single-peaked eccentricity distribution with mean of ~ 0.3) linked with 
Saturn's migration speed ~ 0.8 AU My^^. 

Subject headings: celestial mechanics — minor planets, asteroids: general — 
planets and satellites: dynamical evolution and stability — solar system: general 



1. Introduction 



The dynamical structure of the Kuiper Belt suggests that the ou ter solar system 



rienced a phase of planetesimal-driven migration in its early histo ry (IFernandez fc Iplll984 



expe- 



Malhotralll993l . Il995l : iHahn fc Malhotralll999l : iLevison et al.ll2008f ). Pluto and other Kuiper 
belt objects that are trapped in mean motion resonances (MMRs) with Neptune are explained 
by the outward migration of Neptune due to interactions with a more massive primordial 
planetesimal disk in the outer regions of the solar system (jMalhotralll993l . Il995l ). In addi- 
tion, the so-called the scattered disk of the Kui per belt can also be explained by the outward 
migration of Neptune (IHahn fc Malhotrall2005l ). or by the effects of a hi gh eccentricity phas e 
of ice giant planet evolution during the outward migration of Neptune (ILevison et al.ll2008l ). 
The basic premise of planetesimal-driven migration is that the giant planets formed in a 
more compact configuration than we find them today, and that they were surrounded by 
a massive (~ 50 M®) disk of u naccreted icy planetesim als that was the progenitor of the 
currently observed Kuiper belt (IHahn fc Malhotra]ll999l ). When planetesimals are preferen- 
tially scattered either inward (toward the Sun) or outward (away from the Sun), net orbital 
angular momentum is transfe rred between the disk and the large body , causing a drift in the 
large body's semimajor axis (IFernandez fc Iplll984l : iKirsh et al.ll2009l ). In many simulations 
of giant planet migration, icy planetesimals are preferentially scattered inward by each of the 
three outer giant planets (Saturn, Uranus, and Neptune) causing these planets to migrate 
outward. Due to Jupiter's large mass, planetesimals that encounter Jupiter are preferentially 
ejected out of the solar system, leading to a net loss of mass from the solar system and an 
inward migration of Jupiter. 



Planetesimal-driven giant planet migration has been suggested as a cause of the Late 
Heavy Bombardment (LHB) (iGomes et al.l l2005l : IStrom et al.ll2005l ). however the link be- 
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tween these two events has yet t o be definitively estabhshed f IChapman et aLll2007l : ICuk et al. 
2OIOI : iMalhotra fc Stronj I2OIOI ) . Such migration would have enhanced the impact fiux of 
both asteroids and comets onto the terrestrial planets in two ways. First, many of the icy 
planetesimals scattered by the giant planets would have crossed the orbits of the terrestrial 
planets. Second, as the giant planets migrated, locations of mean motion and secular res- 
onances would have swept across the asteroid belt, raising the eccentricities of asteroids to 
planet-crossing values. 



Recently, iMinton &: Malhotral (120091 ) showed that the patterns of depletion observed in 
the asteroid belt are consistent with the effects of sweeping of resonances during the migration 
of the giant planets. The Jupiter-facing sides of some of the Kirkwood gaps (regions of 
the asteroid belt that are nearly empty due to strong jovian mean motion resonances) are 
depleted relative to the Sun-facing sides, as would be expected due to the inward migration 
of Jupiter and the associated inward sweeping of the jovian mean motion resonances. The 
region within the inner asteroid belt between semimajor axis range 2.1-2.5 AU also has 
excess depletion relative to a model asteroid belt that was uniformly populated and then 
subsequently sculpted by the gravitational perturbations of the planets over 4 Gy, as would 
be expected due to the outward migration of Saturn and the associated inward sweeping of a 
strong secular resonance, the so-called uq resonance, as explained below. In our 2009 study, 
we concluded that the semimajor axis distribution of asteroids in the main belt is consistent 
with the inward migration of Jupiter and outward migration of Saturn by a mounts propose d 



in previous studies based on the Kuiper belt resonance structure(e.g., iMalhotral Il995l ). 
However, in that study the migration timescale was not strongly constrained, because only 
the relative depletion of asteroids in nearby semimajor axis bins could be determined, not 
their overall level of depletion. In the present paper, we explore in more detail the effect that 
planet migration would have had on the asteroid belt due to asteroid eccentricity excitation 
by the sweeping of the uq secular resonance. From the observed eccentricity distribution of 
main belt asteroids, we find that it is possible to derive constraints on the secular resonance 
sweeping timescale, and hence on the migration timescale. 

Secular resonances play an important role in the evolution of the main asteroid belt. 
The inner edge of the belt nearly coincides with the i/g secular resonance which is defined 
by (7 ~ (76) where g is the rate of precession of the longitude of pericenter, w, of an asteroid 
and qq is the sixth eigenfrequency of the solar system planets (approximately the rate of 
precession of Saturn's longitude of pericenter). The z/g resonance , is important for the de - 
livery of Near Earth Astero ids (NEAs) to the inner solar system ( ISchoU fc Froeschlelll99ll ). 



Williams fc Faulknerl (119811 ) showed that the location o 



surfaces in a — e — sini space, and iMilani fc Knezevid (jl990l ) showed that those surfaces 



the z/g resonance actually forms 



approximately define the "inner edge" of the main asteroid belt. However, as mentioned 
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above, iMinton fc Malhotral ( 120091 ) found that, with the planets in their present configura- 
tion, planetary perturbations over the age of the solar system cannot fully account for the 
detailed orbital distribution of the asteroi ds in the inner astero i d belt . The pattern of excess 
depletion in inner asteroid belt noted by iMinton fc Malhotral ( l2009l ) is consistent with the 
effect of the inward sweeping of the secular resonance. In general, the direction of motion 
of the is anticorrelated with that of Saturn, so an inward sweeping of the would be 
produced by an outwardly migrating Saturn. 

Sweeping, or scanning, secular resonances have been analyzed in a number of previous 
works. Sweeping secular resonances due to the changing quadrupole moment of the Sun 
during solar spin-down have been expl ored as a possible mechanism for explaining the ec- 
centricity and inclination of Mercury (jWard et al.lll976l ). Secular resonance sweeping due 
to the effects of the dissipating solar nebula just after planet formation has also been in- 
vestigated as a possible mechanisni for excitin g the orbital eccentricities of Mars and of the 
asteroid belt (lHeppenheimerlll980l : IWardlll98ll ) . The dissipating massive gaseous solar nebula 
would have altered the secular frequencies of the solar system planets in a time-dependent 
way, causing locations of secular resonances to possibly sweep across the inner solar system, 
thereby exciting asteroids in to the eccentric and inc lined orbits that are observed today. This 
mechanism was revisited by iNagasawa et al.l (120001) . who in c orpor ated a more sophisticated 
treatment of the nebular dispersal. However, lO'Brien et al.l (120071 ) have argued that the ex- 
citation (and clearing) of the primordial asteroid belt was unlikely due to secular resonance 
sweeping due to the dispersion of the solar nebula. 

The special case of asteroids on initially circul ar orbits being swept by the z/g and uiq 
resonances has been investigated by iGomed (119971 ). In this paper, we consider the more 
general case of non-zero initial eccentricities; our analysis yields qualitatively new results 
and pro vides new in si ghts into t he dyna r nical h istory of the asteroid belt. This extends the 
work of I Ward et al.l (119761 ) and iGomed (119971 ) in developing analytical treatments of the 
effects of sweeping secular resonances on asteroid orbits. In doing so, we have developed an 
explicit relationship between the migration rate of the giant planets, the initial eccentricity 
of the asteroid and its initial longitude of perihelion, and the final eccentricity of the asteroid 
after the passage of the resonance. We show that for initially non-zero asteroid eccentricity, 
the sweeping of the uq resonances can either increase or decrease asteroid eccentricities. 
Examining the orbits of observed main belt asteroids we find evidence for a double-peaked 
eccentricity distribution; this supports the case for a history of uq sweeping. Quantitative 
comparison of our analytical theory with the semimajor axis and eccentricity distribution of 
asteroids yields new constraints on the timescale of planet migration. 



We note that although our analysis is carried out in the specific context of the sweep- 
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ing z/g resonance during the phase of planetesimal- driven migration of Jupiter and Saturn, 
the techniques developed here may be extended to other similar problems, for example, 
the sweeping of the inclination- node z/ig resonance in the main asteroid belt, the secu- 
lar resonance in the Kuiper belt, and farther afield, the sweeping of secular resonances in 
circumstellar or other astrophysical disks. 



2. Analytical theory of a sweeping secular resonance 

We adopt a simplified model in which a test particle (asteroid) is perturbed only by 
a single resonance, the z/g resonance. We use a system of units where the mass is in solar 
masses, the semimajor axis is in units of AU, and the unit of time is (27r)~^y. With this 
choice, the gravitational constant, G, is unity. An asteroid's secular pertur bations close to a 



secular resonance can be described by the following Hamiltonian function ( lMalhotralll998l ): 



Hsec= -gO-J + l^V2J COs{Wp- ZU), (1) 

where Wp = Qpt + /3p describes the phase of the p-th eige nmode of the linearized ec centricity- 



pericenter secular theory for the Solar system planets (IMurray fc Dermottlll999l ). gp is the 
associated eigenfrequency, zu is the asteroid's longitude of perihelion, J = ^/a (l — \/l — e"^) 
is the canonical generalized momentum which is related to the asteroid's orbital semimajor 
axis a and eccentricity e; —w and J are the canonically conjugate pair of variables in this 
1-degree-of-freedom Hamiltonian system. The coefficients Qq and e are given by: 

^0 = Yl "i^3/2 («i ' (2) 



3 



^ = i^E«?^3/2K-K-^f' (3) 

j 

where the subscript j refers to a planet, is the amplitude of the gp mode in the j*'^ 
planet's orbit, aj = mm{a/aj,aj/a}, rrij is the ratio of the mass of planet j to the Sun, 
and &3/2(Q!i) and b^^i^^j) ^''^^ Laplace coefficients; the sum is over all major planets. The 
summations in equations ([2])-(j3]) are over the 8 major planets, for the greatest accuracy; 
however, we will adopt the simpler two-planet model of the Sun- Jupiter-Saturn in ^ in 
which case we sum over only the indices referring to Jupiter and Saturn; then Qp is an 
eigenfrequency of the secular equations of the two planet system. 

With fixed values of the planetary masses and semimajor axes, go, gp and e are constant 
parameters in the Hamiltonian given in equation ([1]). However, during the epoch of giant 
planet migration, the planets' semimajor axes change secularly with time, so that go, gp and 
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e become time-dependent parameters. In the analysis below, we neglect the time-dependence 
of and e, and adopt a simple prescription for the time-dependence of (see equation [8] 
below). This approximation is physically motivated: the fractional variation of Qq and e for 
an individual asteroid is small compared to the effects of the "small divisor" qq — during 
the 1^6 resonance sweeping event. 

It is useful to make a canonical transformation to new variables (0, P) defined by the 
following generating function, 

J^(-ro, P, t) = {wp{t) - w)P (4) 

Thus, = dT jdP = {wp(t) — w) and J = —dT/dw = P. The new Hamiltonian function is 
Hsec = Hsec + dJ^ / dt, 

Hsec = {wp{t) - go)J + eV2J cos0, (5) 

where we have retained J to denote the canonical momentum, since P = J. It is useful to 
make a second canonical transformation to canonical eccentric variables, 



x = -\/2Jcos0, y = — "\/2Jsin 0, (6) 

where x is the canonical coordinate and y is the canonically conjugate momentum. The 
Hamiltonian expressed in these variables is 

2;^ -|- y'^ 

Hsec = {wp{t) - go) — h ex. (7) 

As discussed above, during planetary migration, the secular frequency gp is a slowly 
varying function of time. We approximate its rate of change, gfp = 2A, as a constant, so that 

Wpi^) = gp,o + 2At. (8) 



We d efine t = as the epoch of exact resonance crossing, so that gp^ = go (cf. I Ward et al. 



19761 ). Then, Wp{t) — go = 2At, and the equations of motion from the Hamiltonian of 
equation ([7]) can be written as: 

X = 2Xty, (9) 

y = -2Xtx-e. (10) 

These equations of motion form a system of linear, nonhomogenous differential equations, 
whose solution is a linear combination of a homogeneous and a particular solution. The 
homogeneous solution can be found by inspection, giving: 

a;fe(t) = ci cos At^ + C2 sin At^, (11) 
Vhit) = — Ci sin At^ + C2 cos Xt^, (12) 
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where ci and C2 are constant coefficients. We use the method of variation of parameters to find 
the particular solution. Accordingly, we replace the constants Ci and C2 in the homogeneous 
solution with functions A{t) and B{t), to seek the particular solution of the form 



Xp{t) = A{t) cos + B{t) sin 
yp{t) = -A{t) sin Xt^ + B{t) cos At . 

Substituting this into the equations of motion we now have: 

icosAt^ + EsinAt^ = 0, 



therefore 



-A sin At + B cos At^ 



A = esin At^, 
B = -ecos At^. 



(13) 
(14) 



(15) 
(16) 



(17) 
(18) 



Equations f lT7|l and f ITH]) do not have a si mple closed-form solution, but their solution can be 
expressed in terms of Fresnel integrals ( IZwillingerl 119961 ) . The Fresnel integrals are defined 
as follows: 



S{t) 

m 



sin t'^dt' 



cos t'^rft'. 



and have the following properties: 



S{-t) = -S{t), 
Ci-t) = -C(t), 

S{oo) = C{oo) 



Therefore 



Ait) ^ -J=S {tVlXl 

Bit) = --^C {ty^\ 



(19) 
(20) 



(21) 
(22) 

(23) 



(24) 
(25) 
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We denote initial conditions with a subscript i, and write the solution to equations (Q 
and ffTOj) as 

x{t) = Xi cos [A {f -t^i)]+ Vi sin [A {f - t-)] 

+ ^= [{S - S,) cos Xt^ -{C- Q) sin Xf] , (26) 

yit) = -X, sin [A - t^)] + y, cos [A (t^ - ^2)] 

[{C - d) cos Xt^ + {S- S,) sin At^] . (27) 



Because the asteroid is swept over by the secular resonance at time t = 0, we can calculate 
the changes in x, y by letting ti = —tj and evaluating the coefficients Cj, C/, Si, Sj far from 
resonance passage, i.e., for t/-\/[Af ^ 1, by use of equation fl23|l . Thus we find 



Xf = Xi + e^l -—- [cosXt'^ — sinXt'^] , (28) 



yj = y,-ej^ [cos At^ + sin At?] . (29) 



The final value of J long after resonance passage is therefore given by 

1 71 s'^ I TT 

= o (^i + + ^ + ^a/ T5ru[^i(cos At,2 - sin Xtf) 



^ ■ ■ 2|A| Y 2|A| 

yi{cos Xtl + sin At?)] 



Tre^ 211 Ji , , ,9 TT, 

= + ^ + ^ cos(0, - At? - -). (30) 

With a judicious choice of the initial time, tj, and without loss of generality, the cosine in 
the last term becomes cosw,-, and therefore 



2p^ V^^^^^^'' 

The asteroid's semimajor axis a is unchanged by the secular perturbations; thus, the 
changes in J reflect changes in the asteroid's eccentricity e. For asteroids with non-zero initial 
eccentricity, the phase dependence in equation (PT]) means that secular resonance sweeping 
can potentially both excite and damp orbital eccentricities. We also note that the magnitude 
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of eccentricity change is inversely related to the speed of planet migration. In linear secular 
theory (equation ([1])), eccentricity and inclination are decoupled, an d therefore the effect 



of the sweeping z/g does not depend on the inclination. However, as IWilliams fc Faulkner 



(Il98ll ) showed, the location of the z/g does depend on inclination, but the dependence is weak 
for typical inclinations of main belt objects. Nevertheless, there are populations of main belt 
asteroids at high inclination (such as the Hungaria and Phocaea families), and an analysis of 
secular resonance sweeping that incorporates coupling between eccentricity and inclination 
would be valuable for understanding the effects of planet migration on these populations; we 
leave this to a future investigation. 

For small e, we can use the approximation J ~ |-y/ae^. Considering all possible values 
of cosGJj G {—1, +1}, an asteroid with initial eccentricity Cj that is swept by the resonance 
will have a final eccentricity in the range Cmm to Cmax, where 



and 

I 'IT 

(33) 



TT 



Equations (l3T]) - (l33l) have the following implications: 



1. Initially circular orbits become eccentric, with a final eccentricity Se- 

2. An ensemble of orbits with the same a and initial non-zero e but uniform random 
orientations of pericenter are transformed into an ensemble that has eccentricities in 
the range emin to Cmax] this range is not uniformly distributed because of the cost^j 
dependence in equation f l3T]) . rather the distribution peaks at the extreme values (see 
Fig. [1] below). 

3. An ensemble of asteroids having an initial distribution of eccentricities which is a single- 
peaked Gaussian (and random orientations of pericenter) would be transformed into 
one with a double-peaked eccentricity distribution. 



3. Application to the Main Asteroid Belt 

In light of the above calculations, it is possible to conclude that if the asteroid belt were 
initially dynamically cold, that is asteroids were on nearly circular orbits prior to secular 
resonance sweeping, then the asteroids would be nearly uniformly excited to a narrow range of 
final eccentricities, the value of which would be determined by the rate of resonance sweeping. 
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Because asteroids having eccentricities above planet-crossing values would be unlikely to 
survive to the present day, it follows that an initially cold asteroid belt which is uniformly 
excited by the uq sweeping will either lose all its asteroids or none. On the other hand, an 
initially excited asteroid belt, that is a belt with asteroids that had non-zero eccentricities 
prior to the sweeping, would have asteroids' final eccentricities bounded by equation (132|) . 
allowing for partial depletion and also broadening of its eccentricity distribution. In this 
section, we apply our theoretical analysis to the problem of the resonance sweeping through 
the asteroid belt, and compare the theoretical predictions with the observed eccentricity 
distribution of asteroids. 



3.1. Parameters 



In order to apply the theory, we must find the location of the resonance as a function 
of the semimajor axes of the giant planets orbits, and also obtain values for the parameter e 
(equation [3]), for asteroids with semimajor axis values in the main asteroid belt. The location 
of the resonance is defined as the semimajor axis, a^^g, where the rate, go (equation [2]) , of 
pericenter precession of a massless particle (or asteroid) is equal to the eigenfrequency of 
the solar system. In the current solar system, the ge frequency is associated with the secular 
mode with the most power in Saturn's eccentricity-pericenter variations. During the epoch 
of planetesimal-driven planet migratio n, Jupiter migrated by only a small amount but Sat 



urn l ikely migrated significantly more ( iFernandez fc Iplll984t lMalhotralll995l : iTsiganis et al. 



20051 ). so we expect that the variation in location of the uq secular resonance is most sensitive 
to Saturn's semimajor axis. We therefore adopt a simple model of planet migration in which 
Jupiter is fixed at 5.2 AU and only Saturn migrates. We neglect the effects of the ice giants 
Uranus and Neptune, as well as secular effects due to the previously more massive Kuiper 
belt and asteroid belt. With these simplifications, the g^ frequency varies with time as Sat- 
urn migrates, so we parametrize g^ as a function of Saturn's semimajor axis. In contrast 
with the variation of gg, there is little variation of the asteroid's apsidal precession rate, go, 
as Saturn migrates. Thus, finding a^^ is reduced to calculating the dependence of gQ on 
Saturn's semimajor axis. 

To calculate g^ as a function of Saturn's semimajor axis, we proceed as follows. For 
fixed planetary semimajor axes, the Laplace-Lagrange secular theory provides the secular 
frequencies and orbital element variations of the planets. This is a linear perturbation theory, 
in which the disturbing function is truncated to s ecular terms of second order in eccentricity 
and first order in mass (IMurray fc Dermottlll999l ). In the planar two-planet case, the secular 
perturbations of planet j, where j = 5 is Jupiter and j = 6 is Saturn, are described by the 
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following disturbing function: 



Ajje^j + AjkC^ee cos(g75 - zuq 



where n is the mean motion, and A is a matrix with elements 



A 



for j 



5,6, k = 6, 5, and j ^ k; ajk = min{aj/afc, ak/aj}, and 



a 



jk 



1 • % > (^k 



(34) 



(35) 
(36) 



(37) 



The secular motion of the planets is then described by a set of linear differential equations 
for the eccentricity vectors, ej(sinccTj, costu^) = {hj,kj). 



hi 



p=5 



p=5 



(38) 



For fixed planetary semimajor axes, the coefficients are constants, and the solution is given 
by a linear superposition of eigenmodes: 



{hj, kj} = ^ Ej^^{cos(5(pt + Pp), sm{gpt + /3p)}, 



(39) 



where gp are the eigenfrequencies of the matrix A and Ej^^ are the corresponding eigen- 
vectors; the amplitudes of the eigenvectors and the phases (3p are determined by initial 
conditions. In our 2-planet model, the secular frequencies and g^ depend on the masses 
of Jupiter, Saturn, and the Sun and on the semimajor axes of Jupiter and Saturn. 

For the current semimajor axes and eccentricities of Jupiter and Saturn the Laplace- 
Lagrance theory gives frequency values = 3.7" y~^ and = 2 2 .3^^ y~ ^, which are lower 
than the mo re accurate va l ues given bv iBrouwer &: van WoerkomI (jl950l ) by 14% and 20%, 
respectively (jLaskarlll988l ). iBrouwer fc van Woerkom Jl950h achieved their more accurate 
solution by incorporating higher order terms in the disturbing function involving 2A5 — 5Xq, 
which arise due to Jupiter and Saturn's proximity to the 5:2 resonance (the so-called "Great 
Inequality"). By doing an accurate numerical analysis (described below), we found that the 
effect of the 5:2 resonance is only important over a very narrow range in Saturn's semimajor 
axis. More significant is the perturbation owing to the 2:1 near-resonance of Jupiter and 



-VI- 



Saturn. iMalhotra et al.l (jl989[ ) developed corrections to the Laplace-Lagrange theory to 
account for the perturbations from n^\: n resonances in the context of the Uranian sateUite 
system. Applying that approach to our problem, we find that the 2:1 near-resonance between 
Jupiter and Saturn leads to zeroth order corrections to the elements of the A matriji0. 
Including these corrections, we determined the secular frequencies for a range of values of 
Saturn's semimajor axis; the result for is shown in Fig. (dashed line). 

We also calculated values for the eccentricity-pericenter eigenfrequencies by direct nu- 
merical integration of the full equations of motion for the two-planet, planar solar system. 
In these integrations, Jupiter's initial semimajor axis was 5.2 AU, Saturn's semimajor axis, 
ag, was one of 233 values in the range 7.3-10.45 AU, initial eccentricities of Jupiter and 
Saturn were 0.05, and initial inclinations were zero. The initial longitude of pericenter and 
mean anomalies of Jupiter were w^^i = 15° and X^^i = 92°, and Saturn were cug^i = 338° and 
A6,j = 62.5°. In each case, the planets' orbits were integrated for 100 myr, and a Fourier 
transform of the time series of the {hj, kj} yields their spectrum of secular frequencies. For 
regular (non-chaotic) orbits, the spectral frequencies are well defined and are readily identi- 
fied with the frequencies of the secular solution. The frequency as a function of Saturn's 
semimajor axis was obtained by this numerical analysis; the result is shown by the solid line 
in Fig. Hi. 

The comparison between the numerical analysis and the Laplace-Lagrange secular the- 
ory indicates that the linear secular theory, including the corrections due to the 2:1 near- 
resonance, is an adequate model for the variation in g^ as a function of ag. We adopted 
the latter for its convenience in the needed computations. The value of function of 

Saturn's semimajor axis was thus found by solving for the value of asteroid semimajor axis 
where go = g^; go was calculated using equation ([2]) and g^ is the eigenfrequency associated 
with the p = 6 eigenmode (at each value of Saturn's semimajor axis). The result is shown 
in Fig. [2)3. 

We also used the analytical secular theory to calculate the eigenvector components E^^^ 
in the secular solution of the 2-planet system, for each value of Saturn's semimajor axis. 
We adopted the same values for the initial conditions of Jupiter and Saturn as in the direct 
numerical integrations discussed above. Finally, we computed the values of the parameter 
e at each location Oj^g of the secular resonance. The result is plotted in Fig. [31 Despite the 



^Corrections due to near-resonances are of order e^""^-*, where e is eccentricity and n is the order of the 
resonance. The 5:2 is a third order resonance, so its effect is O(e^). The 2:1 is a first order resonance, so that 
its effect does not depend on e. Therefore the discrep ancy between hnear theory and numerical analysis (or 
the higher order theory of iBrouwer fc van Woerkoml ) arising from the Great Inequality would be much less 
if Jupiter and Saturn were on more circular orbits, but the effect due to the 2:1 resonance would remain. 
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complexity of the computation, the result shown in Fig. |3] is approximated well by a simple 
exponential curve, e ~ 3.5x 10~^ exp(2at,g/AU), in the semimajor axis range 2 < a^,g/AU < 4. 



3.2. Four test cases 

We checked the results of our analytical model against four full numerical simulations of 
the restricted four-body problem (the Sun- Jupiter-Saturn system with test particle asteroids) 
in which the test particles in the asteroid belt are subjected to the effects of a migrating 
Saturn. The numerical integratioii was performed with an implementation of a symplectic 



mapping (jWisdom fc Holmanlll99ll : ISaha fc Tremainelll992l ). and the integration stepsize was 



0.01 y. Jupiter and Saturn were the only massive planets that were integrated, and their 
mutual gravitational influence was included. The asteroids were approximated as massless 
test particles. The current solar system values of the eccentricity of Jupiter and Saturn 
were adopted and all inclinations (planets and test particles) were set to zero. An external 
acceleration was applied to Saturn to cause it to migrate outward linearly and smoothly 
starting at 8.5 AU at the desired rate. As shown in Fig. [2]3, the resonance was located at 
3 AU at the beginning of the simulation, and swept inward past the current location of the 
inner asteroid belt in all simulations. In each of the four simulations, 30 test particles were 
placed at 2.3 AU and given different initial longitudes of pericenter spaced 12° apart. The 
semimajor axis value of 2.3 AU was chosen because it is far away from the complications 
arising due to strong mean motion resonances. The only parameters varied between each of 
the four simulations were the initial osculating eccentricities of the test particles, e^, and the 
migration speed of Saturn, dg- The parameters explored were: 

a) d = 0.2, de = 1.0 AU My"^; 

b) ei = 0.2, de = 0.5 AU My'^ 

c) ei = 0.1, de = 1.0 AU My"^; 

d) ei = 0.3, de = 1.0 AU My"\ 

These values were chosen to illustrate the most relevant qualitative features. The migration 
rates of 0.5 AU My~^ and 1.0 AU My~^ are slow enough so that the change in eccentricity is 
substantial, but not so slow that the non-linear effects at high eccentricity swamp the results. 
These test cases illustrate both how well the analytical model matches the numerical results, 
and where the analytical model breaks down. 
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Two aspects of the analytical model were checked. First, the perturbative equations 
of motion, equations iQ and (fTOl) . were numerically integrated, and their numerical solu- 
tion compared with the numerical solution from the direct numerical integration of the full 
equations of motion. For the perturbative solution, we adopted values for A that were ap- 
proximately equivalent to the values of in the full numerical integrations. Second, the 
eccentricity bounds predicted by the analytical theory, equation (l32l) . were compared with 
both numerical solutions. The results of these comparisons for the four test cases are shown 
in Fig. m We find that the analytically predicted values of the maximum and minimum final 
eccentricities (shown as horizontal dashed lines) are in excellent agreement with the final 
values of the eccentricities found in the numerical solution of the perturbative equations, 
and in fairly good agreement with those found in the full numerical solution. Not surpris- 
ingly, we find that the test particles in the full numerical integrations exhibit somewhat more 
complicated behavior than the perturbative approximation, and equation fl5^ somewhat un- 
derpredicts the maximum final eccentricity: this may be due to higher order terms in the 
disturbing function that have been neglected in the perturbative analysis and which become 
more important at high eccentricity; effects due to close encounters with Jupiter also become 
important at the high eccentricities. 



3.3. Comparison with observed asteroid eccentricities 



Does the eccentricity distribution of main belt asteroids retain features corresponding to 
the effects of the uq resonance sweeping? To answer this question, we need to know the main 
belt eccentricity distribution free of observational bias, and also relatively free of the effects of 
~ 4 Gy of collisional evolution subsequent to the effects of planetary migration. We therefore 
obtained the proper elements of the observationally complete sa mple of asteroids with ab - 



solute magnitude H < 10.8 from the AstDys online data service jKnezevic fc Milanilboosh: 



we ex cluded from this set the members of c ollisional families as identif ied by iNesvorny et al. 



(120061 ). These same criteria were adopted in lMinton Sz Malhotral (120101 ) in a study of the long 
term dynamical evolution of large asteroids. This sample of 931 main belt asteroids is a good 
approximation to a complete set of large asteroids that have been least perturbed by either 
dynamical evolution or collisional evolution since the epoch of the last major dynamical event 
that occurred in this region of the solar system; therefore this sample likely preserves best the 
post-migration orbital distribution of the asteroid belt. The proper eccentricity distribution 
of these asteroids is shown in Fig. O This distribution has usually been described in the lit 



erature by simply quoting its mean value (and sometimes a dispersion) (IMurray fc Dermott 



19991 : lO'Brien et al.ll2007l ). Our best fit single Gaussian distribution to this data has a mean, 
fie and standard deviation, (Tg, given by /ie = 0.135 ±0.00013 and Ue = 0.0716 ±0.00022, and 
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is plotted in Fig. |5l However, we also note (by eye) a possible indication of a double-peak 
feature in the observed population. Our best fit double Gaussian distribution (modelled as 
two symmetrical Gaussians with the same standard deviation, but different mean values) to 
the same data has the following parameters: 

= 0.0846 ±0.00011, 
/i^ s = 0.185 ±0.00012, 
a'^ = 0.0411 ± 0.00020. 

More details of how these fits were obtained are described in Appendix \^ where we also 
discuss goodness-of-fit of the single and double gaussian distributions. A Kolmogorov- Sinai 
test shows that there is only a 4.5% probability that the observed eccentricities actually are 
consistent with the best-fit single Gaussian distribution, and a 73% probability that they are 
consistent with the best-fit double Gaussian distribution; while the double gaussian is appar- 
ently a better fit to the data compared to the single gaussian, a K-S probability of only 73% 
is quite far from a statistically significa nt level of confidence. A dip test for mult i- modality 



in the observed eccentricity distribution (iHartigan fc Hartiganlll985l ) is also inconclusive; the 



likely reason for this is that our data sample is not very large, and the separation between 
the two putative peaks is too small in relation to the dispersion. Nevertheless, we bravely 
proceed in the next section with the implications of a double-Gaussian eccentricity distribu- 
tion, with the caveat that some of the conclusions we reach are based on this statistically 
marginal result. 



4. A constraint on Saturn's migration rate 

By relating the Qq secular frequency to the semimajor axis of Saturn, can be re- 
lated to the migration rate of Saturn, dg- In this section only the effects of the sweeping 
z/g secular resonance will be considered, and effects due to sweeping jovian mean motion 
resonances will be ignored. Because Jupiter is thought to have migrated inward a much 
smaller distance than Saturn migrated outward during planetesimal-driven migration, the 
effects due to migrating jovian MMR s were likely confined to narrow regions near strong 



resonances (IMinton fc Malhotrall2009l ). In the inner asteroid belt between 2.1-2.8 AU, these 
would include the 3:1 and 5:2 resonances, currently located at approximately 2.5 AU and 
2.7 AU, respectively. As shown in Fig. |2]3, plausible parameters for the outward migration 
of Saturn would have allowed the z/g to sweep across the entire inner asteroid belt. Therefore 
the z/g resonance would have been the major excitation mechanism across the 2.1-2.8 AU 
region of the main belt (and possibly across the entire main asteroid belt, depending on 
Saturn's pre-migration semimajor axis) during giant planet migration. 
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We used the results of our analytical model to set limits on the rate of migration of 
Saturn assuming a linear migration profile, with the caveat that many important effects 
are ignored, such as asteroid- Jupiter mean motion resonances, and Jupiter-Saturn mean 
motion resonances (with the exception of the 2:1 resonance). We have confined our analysis 
to only the region of the main belt spanning 2.1-2.8 AU. Beyond 2.8 AU strong jovian 
mean motion resonance become more numerous. Due to the high probability that the icy 
planetesimals driving planet migration would be ejected from the solar system by Jupiter, 
Jupiter likely migrated inward. The migration of Jupiter would have caused strong jovian 
mean motion resonances to sweep the asteroid belt, ca.using additional depletion beyond 



that of the sweeping uq resonance (IMinton fc Malhotral l2009l ). A further complication is 



that sweeping jovian mean motion resonances may have also trapped icy plan etesimals that 



enter ed the asteroid belt region from their source region beyond Neptune f lLevison et al. 



20091 ). The effects of these complications are reduced when we consider only the inner 
asteroid belt. From Fig. |2]d, we find that the would have swept the inner asteroid belt 
region between 2.1-2.8 AU when Saturn was between ~ 8.5-9.2 AU. Therefore the limits 
on de that we set using the inner asteroid belt as a constraint are only applicable for this 
particular portion of Saturn's migration history. 

Our theoretically estimated final eccentricity as a function of initial asteroid semimajor 
axis and eccentricity is shown in Fig. [6] for three different adopted migration rates of Saturn. 
The larger the initial asteroid eccentricities, the wider the bounds in their final eccentricities. 
If we adopt the reasonable criterion that an asteroid is lost from the main belt when it 
achieves a planet-crossing orbit (that is, crossing the orbits of either Jupiter or Mars) and 
that initial asteroid eccentricities were therefore confined to < 0.4, then from Fig. |6] Saturn's 
migration rate must have been clq > 0.15 AU My~^ when the uq resonance was sweeping 
through the inner asteroid belt. Our results indicate that if Saturn's migration rate had 
been slower than 0.15 AU My~^ when it was migrating across ~ 8.5-9.2 AU, then the inner 
asteroid belt would have been completely swept clear of asteroids by the uq resonance. 

In light of our analysis and the observed dispersion of eccentricities in the asteroid belt 
(Fig. |5]), we can also immediately conclude that the pre- migration asteroid belt between 
2.1 and 2.8 AU had significantly non-zero eccentricities. This is because, as discussed at 
the top of §31 an initially cold asteroid belt swept by the z/g resonance would either lose 
all its asteroids or none, and very low initial eccentricities would result in final asteroid 
eccentricities in a very narrow range of values (cf. equation [32]) . in contradiction with the 
fairly wide eccentricity dispersion that is observed. This conclusion supports recent results 
from studies of planetesimal accretion and asteroid and pl anet formation th at the asteroids 



were modestly excited at the end of their formation (e.g., iPetit et al.ll2002[ ) 
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In Appendix |X] we show that the double-gaussian distribution is a shghtly better fit 
to the main belt asteroid eccentricity distribution, but the statistical tests do not rule out 
a single-peaked distribution. We boldly proceed with considering the implications of the 
double-peaked eccentricity distribution to further constrain the migration rate of Saturn, 
with the caveat that these results can only be said to be consistent with the observations, 
rather than uniquely constrained by them. 

If the pre-sweeping asteroid belt had a Gaussian eccentricity distribution, then the lower 
peak of the post-sweeping asteroid belt should be equal to the lower bound of equation (15^ . 
We use the analytical theory to make a rough estimate of the parameter A (and hence de) 
that would yield a final distribution with lower peak near 0.09 and upper peak near 0.19 
(which is similar to the best-fit double Gaussian in Fig. . 



Applying equation fl5^ . we see that there are two possible solutions: (cj) = 0.14,6^ = 
0.05 and (cj) = 0.05,6e = 0.14. A corresponding migration rate of Saturn can be esti- 
mated from the value of 6e using equation (!33|) . and the parameter relationships plotted 
in Figs. |2] and [31 The former solution (5e = 0.05) requires a migration rate for Saturn of 
dg = 30 AU My~^. We mention this implausible solution here for completeness, but we 
will not discuss it any further. The latter solution {6e = 0.14) requires a migration rate for 
Saturn of dg = 4 AU My~^. We dub this solution the "cold belt" solution. This rate is com- 
pa rable to the rates of planet migration found in the "Jumping Jupiter" scenario proposed 



by lBrasser et al.l (j2009[ ). A third solution exists if we consider that eccentricities in the main 
belt are restricted by the orbits of Mars and Jupiter on either side, such that stable asteroid 
orbits do not cross the planetary orbits. This hmits asteroid eccentricities to values such that 
neither the aphelion of the asteroid crosses the perihelion distance of Jupiter, nor the peri- 
helion of the asteroid crosses the aphelion distance of Mars. Maximum asteroid eccentricity 
is therefore a function of semimajor axis, where e-max = niin(l — Quars/O' — 1, <iJupiter/0' — 1), 
where Q and q are planet aphelion and perihelion, respectively, and a is the semimajor axis 
of the asteroid. In this case, an initial single Gaussian eccentricity distribution with a mean 
greater than ~ 0.3 would be severely truncated, therefore we need only fit the lower peak of 
the double Gaussian distribution at e = 0.09. Applying equation fl32|) . we find that 5^ = 0.21 
provides a good fit. The corresponding migration rate of Saturn is dg = 0.8 AU My~^. We 
dub this solution the "hot belt" solution. 

We illustrate the two possible solutions for an ensemble of hypothetical asteroids having 
semimajor axes uniformly distributed randomly in the range 2.1 AU to 2.8 AU. In Fig. [7^ the 
initial eccentricity distribution is modeled as a Gaussian distribution with a mean (cj) = 0.05 
and a standard deviation of 0.01. This initial standard deviation was chosen so that the final 
standard deviation would be the same as that of the observed main belt. Fig. [7)d shows the 
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eccentricity distribution after z/g resonance sweeping has occurred due to the migration of 
Saturn at a rate of 4 AU My^^. The final distribution was calculated with equation f l3T]) : we 
used values of e shown in Fig. El and the value of A was calculated with the aid of Fig. |5]which 
relates the value of gp to the semimajor axis of Saturn. As expected, when an ensemble of 
asteroids with a single-peaked eccentricity distribution is subjected to the sweeping secular 
resonance, the result is a double-peaked eccentricity distribution. Because of the slight bias 
towards the upper limit of the eccentricity excitation band, proportionally more asteroids 
are found in the upper peak. 

In Fig. [7t the initial eccentricity distribution is modeled as a truncated Gaussian: a 
Gaussian with mean (cj) = 0.4 and standard deviation 0.1, but truncated at the semimajor 
axis-dependent Mars-crossing value. We used equation fl^ to calculate the eccentricity dis- 
tribution of this hypothetical ensemble after resonance sweeping with dg = 0.8 AU My~^. 
Again, allowing that only those asteroids whose final eccentricities are below the Mars- 
crossing value will remain, the resulting post-migration eccentricity distribution is shown in 
Fig. [Ttl. In this case, we find the lower peak at the same eccentricity value as the lower peak 
in the observed main belt distribution (see Fig. 

In both cases of possible solutions (initially cold main belt with (cj) = 0.05 and dg = 
4 AU My"^; and initially hot main belt with (e^) = 0.4 and dg = 0.8 AU My~^), the 
theoretical models yield an excess of asteroids wi th eccentricities grea t er tha n 0.2 than in 



the observed main belt. However, as shown by iMinton fc Malhotral (120101 ). on gigayear 
timescales, the e > 0.2 population of the asteroid belt is dynamically more unstable than 
the e < 0.2 population. Thus, both solutions may be consistent with the observations, as 
post-sweeping dynamical erosion could result in a final eccentricity distribution resembling 
more closely the observed distribution. 

The estimates of Saturn's migration rate quoted above depend strongly on the eccentric- 
ities of the giant planets during their migration. In deriving the above estimates, we adopted 
the present values of the giant planets' orbital eccentricities. The z/g resonance strength coef- 
ficient e (equation [3]) is proportional to the amplitude of the p = 6 mode, which is related to 
the eccentricities of the giant planets (namely Saturn and Jupiter). From equation fl39p . and 
the definition ej(sin tz^j, cos ccj) = {hj,kj), the value of Ej^^ is a linear combination of the 
eccentricities of the giant planets. Because Saturn is the planet with the largest amplitude 
of the p = 6 mode, from equation (I3T!) the relationship between the sweep rate and the 
value of Saturn's eccentricity is approximately Xmin c< eg. Therefore, to increase the limiting 
timescale by a factor of ten would only require that the giant planets' eccentricities were 
~ 0.3x their current value (i.e., e^^Q ~ 0.015). 
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5. Conclusion and Discussion 

Based on the existence of the inner asteroid belt, we conclude that Saturn's migration 
rate must have been > 0.15 AU My~^ as Saturn migrated from 8.5 to 9.2 AU (as the z/g 
resonance migrated from 2.8 AU to 2.1 AU). Migration rates lower than ~ 0.15 AU My~^ 
would be inconsistent with the survival of any asteroids in the inner main belt, as the uq 
secular resonance would have excited asteroid eccentricities to planet-crossing values. This 
lower limit for the migration rate of Saturn assumes that Jupiter and Saturn had their 
current orbital eccentricities; the migration rate limit is inversely proportional to the square 
of the amplitude of the secular mode; if Jupiter and Saturn's eccentricities were ~ 0.3 x 
their current value (i.e., 65^5 ~ 0.015) during the planet migration epoch, the limit on the 
migration rate decreases by a factor of ~ 10. (This caveat also applies to the migration rate 
limits quoted below.) 

Our analysis of secular resonance sweeping predicts that a single-peaked eccentricity dis- 
tribution will be transformed into a double-peaked eccentricity by secular resonance sweep- 
ing. We find that the observed eccentricity distribution of asteroids may be consistent with 
a double-peaked distribution function, although we acknowledge that the statistics are poor. 
We used a double-Gaussian function to model the observed asteroid eccentricity distribution 
to set even tighter, albeit model-dependent, constraints on the migration rate of Saturn. 

We identified two possible migration rates that depend on the pre-migration dynamical 
state of the main asteroid belt. The first, the "cold primordial belt" solution, has an asteroid 
belt with an initial eccentricity distribution modeled as a Gaussian with (cj) = 0.05; Saturn's 
migration rate of 4 AU My~^ yields a final eccentricity distribution consistent with the 
observed asteroid belt. The second, the "hot primordial belt" solution, has an asteroid belt 
with an initial eccentricity distribution modeled as a Gaussian with (ej) = 0.4, cxe = 0.1, but 
truncated above the Mars-crossing value of eccentricity; in this case, Saturn's migration rate 
of 0.8 AU My~^ is generally consistent with the observed asteroid belt. 

Each of these solutions has very different implications for the primordial excitation 
and depletion of the main asteroid belt. The cold belt solution, with dg = 4 AU My~^, 
would lead to little depletion of the asteroid belt during giant planet migration, as the 
resonance would be unable to raise eccentricities to Mars-crossing values. This implies an 
initial dynamically quite cold asteroid belt with not more than ~ 2 times the mass of the 
current main belt; the l atter estimate comes from accounting for dynamical erosion over the 



age of the solar system (IMinton fc Malhotra]l2010l ). 



The hot belt solution, with dg = 0.8 AU My ^, would lead to loss of asteroids directly 
due to excitation of eccentricities above planet-crossing values, by about a factor of ~ 2. In 
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this case, the main asteroid belt was more dynamically excited prior to resonance sweeping 
than we find it today. This implies much greater loss of asteroids prior to z/g-sweeping, as 
the peak of the eccentricity d istribution would be near the Mars-crossing value, and subject 
to strong dynamical erosion (IMinton fc Malhotra]l2010l ). 



Each of these solutions has different implications for the model that the Late Heavy 
Bombardment of the inner solar system is linked to th e epoch of planetesimal-driven giant 
planet migration ( iGomes et al.ll2005l : IStrom et al.ll2005l ). These implications will be explored 
in a future paper. 

We remind the reader that in order to elucidate the effects of uq resonance sweeping, 
we have made a number of simplifying assumptions to arrive at an analytically tractable 
model. These simplifications include neglecting the effects of planets other than Jupiter and 
Saturn, the effects of sweeping jovian mean motion resonances on asteroids, the effects of a 
presumed massive Kuiper belt during the epoch of planet migration, and the self-gravity and 
collisional interactions of a previously more massive asteroid belt. In addition, our analysis 
was carried out in the planar approximation, thereby neglecting any eccentricity-inclination 
coupling effects. These neglected effects can be expected to reduce somewhat the lower limit 
on Saturn's migration speed that we have derived, because in general they would reduce 
the effectiveness of the z/g in exciting asteroid eccentricities. Perhaps more importantly, 
giant planet migration would al so lead to the sweeping of the main asteroid belt by the 
z/i6 inclination secular resonance (jWilliams Sz Faulknerlll98ll ) whose effects could be used to 
infer additional constraints. 



Recently, iMorbidelli et al.l ( l2010l ) found through numerical simulations with slow rates 
of planet migration (e-folding timescales exceeding 5 Myr) that the surviving asteroids in 
the main belt tend to be clumped around mean motion resonances. The semimajor axis 
distribution of survivors is found very different from the observed distribution for asteroids, 
and also that a large proportion of survivors have inclinations above 20°, both inconsistent 
with the current main asteroi d belt; comparisons made w ith the unbiased main belt asteroid 
distributions as described in iMinton Sz Malhotral (|2009[ ). This is likely due to the effects 
of the sweeping uiq inclination-longitude of ascending node secular resonance, analogous to 
the sweeping uq eccentricity-pericenter secular resonance that we analyzed in the present 
study. While the effects of the sweeping uiq resonance are analogous to the z/g, only affecting 
inclinations instead of eccentricities, a full analysis of the asteroid belt inclinations is beyond 
the scope of the present work, but will be explored in a future study. 

A number of other studies have derived limit s on the speed of planetesimal-driven giant 
planet migration. [Murray- Clay fc Chiangl (120051 ) exclude an e— folding migration timescale 
r < 1 My to 99.65% confidence based on the lack of a large observed asymmetry in the 
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population of Kui per belt objects i n the two libration centers of the 2:1 Neptune mean 
motion resonance. iBoue et al.l (j2009[ ) exclude r < 7 My based on the observed obliquity of 



Saturn. The latter lower limit on the migration timescale is slightly incompatible with the 
lower limit on the rate of Saturn's migration of cig > 0.15 AU My~^ we derive based on the 
existence of the inner asteroid belt. One way these can be reconciled is if Saturn's orbital 
eccentricity were a factor ~ 2 smaller than its present value as it migrated from 8.5 AU to 
9.2 AU; then, some mechanism would need to have increased Saturn's eccentricity up to its 
present value by the time Saturn reached its present semimajor axis of ~ 9.6 AU. 
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A. The Main belt eccentricity distribution 

The binned eccentricity distribution may be modeled as a Gaussian probability distri- 
bution function, given by: 
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where a is the standard deviation, /i is the mean, and x is the random variable; in our 
case X is the eccentricity. With an appropriate scaling factor, equation (lAll) can be used 
to model the number of asteroids per eccentricity bin. However, rather than fit the binned 
distribution, we instead perform a least squares fit of the unbinned sample to the Gaussian 
cumulative distribution function given by: 



P{x) 



1 
2 



1 



-erf 



X — II 



(A2) 



2'""V o^f2 

For the eccentricities of our sample of 931 main belt asteroids, the best fit parameters are: 

/ie = 0.135 ±0.00013, 
a, = 0.0716 ± 0.00022. 



We also fit the data to a double- Gaussian distribution 
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The cumulative distribution function for equation (1A3P is 

P2{X) = 1 + 
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For the eccentricities of our sample of 931 main belt asteroids, we performed a least squares 
fit to equation ( ]A4p and obtained the following best-fit parameters: 



e,l 



0.0846 ±0.00011, 



/i^ a = 0.185 ±0.00012, 
aL = 0.0411 ± 0.00020. 



We evaluated the goodness of fit using the Kolmogorov-Smirnov (K-S) test. The K-S 
test determines the probability that two distr ibutions are the s ame, or in our case how well 
our model distributions fit the observed data (jPress et al.lll992l ). The K-S test compares the 
cumulative distribution of the data against the model cumulative distribution function. We 
found that our asteroid sample has a probability of 4.5 x 10~^ that it comes from the best 
fit single Gaussian (equation ( 1A2I) ). but a probability of 0.73 that it comes from the double- 
Gaussian (equation flA4l) ). Therefore, the K-S tests indicate that the double- Gaussian is a 
better fit to the data than the single-Gaussian. 



We performed Hartigan's dip test flHartigan fc Hartiganlll985l ) to test whether the obser- 
vational data is consistent with a multi-peaked distribution. Hartigan's dip test calculates the 
probability that the distribution being tested has a single peak. Applying Hartigan's dip test 
to a given distribution yields in a test statistic; together with the samp le size, the test statis- 
tic is matched to a p- value range in a precomputed table provided by iHartigan fc Hartigan 
(119851 ). The p- value is a measure of the probability that the distribution actually has only 
one peak (the null- hypothesis, for this problem). The smaller the calculated p- value, the 
less likely is the distribution to have a single peak and the more likely it is to have at least 
two significant peaks. A p-value of < 0.05 indicates that the null-hypothesis is very unlikely, 
and that the given distribution has more than one peak. We applied the dip test to the ec- 
centricity distribution of our sample of 931 ma i n bel t asteroids; we determined that the test 
statistic is 0.0107. From lHartigan &: Hartigan! (119851 ). this corresponds to a p-value range of 
0.6 < p < 0.7 (based on a sample size of 1000). This indicates that the Hartigan dip test 
does not rule out the null hypothesis, i.e., a single-peaked distribution cannot be ruled out. 

To further aid the interpretation of this test, we compare this result to that obtained by 
applying Hartigan's dip test to synthetic distributions that were explicitly double-peaked by 
construction. Ten model distributions (of 1000 eccentricity values) were generated from the 
double-gaussian function of equation (1A4I) . with parameter values fii = 0.0846, fi2 = 0.185, 
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and a = 0.0411 (same parameter values as the best-fit for our sample of asteroids). Only the 
random seed was varied between each model distribution. Applying the Hartigan dip test, 
we find that the test statistic ranged between 0.00852 and 0.0126, corresponding to p-values 
between 0.98 and 0.3. This means that, according to the dip test, the null hypothesis - 
i.e., a single-peaked distribution - could not be ruled out for any of the model distributions 
(since none had p < 0.05), despite the fact that they were each drawn from an explicitly 
double-peaked distribution by construction. This illustrates that even with a sample size of 
nearly 1000, the dispersion in eccentricity around the two peaks is too large in comparison 
to the distance between the peaks, that the dip test is not sufficiently sensitive to detect 
the double-peaked underlying distribution. We interpret this to mean that the results of 
both the K-S test and Hartigan's dip test indicate that the main asteroid belt eccentricity 
distribution is consistent with being drawn from a double-peaked distribution, but that this 
cannot be definitely shown. 
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Fig. 1. — The final eccentricity distribution of an ensemble of particles, all having initial 
eccentricity Cj = 0.1 but uniformly distributed values of the phase angle Wi. The effect due 
to the sweeping resonance was modeled using equation fl3T|) . with parameters chosen to 
simulate asteroids at a = 2.3 AU, and with de = 1 AU My~^. 
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Fig. 2. — a) The Qq eigenfrequency as a function of Saturn's semimajor axis, for Jupiter fixed 
at 5.2 AU. The dashed hne shows the result from hnear secular theory, with a correction for 



the e ffect of the near 2:1 mean motion resonance between Jupiter and Saturn (iMalhotra et al. 



19891 ). The solid line shows the result from numerical spectral analysis of 233 solar system 



integrations (see text for explanation). The locations of Jupiter-Saturn MMRs which have 
an effect on the value of Qq are indicated by vertical dotted lines, b) The location of the 
resonance (at zero inclination) as a function of Saturn's orbit. The frequencies and 
Qo were calculated for each value of Saturn's semimajor axis, Og, and then the location a^^ 
was determined by finding where go — 96 = ^- The dashed line shows the result from linear 
secular theory, with a correc tion for the effect of the 2:1 near-MMR between Jupiter and 
Saturn (IMalhotra et al.lll989l ). The solid line was obtained by using the eigenfrequencies 
obtained from spectral analysis of the 233 numerical integrations, as shown in (a). 
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Fig. 3. — The value of the coefficient e defined by equation (jS]) as a function of the zero 
inchnation location of the resonance. The values of E^^^ were calculated using first order 
Laplace-Lagrange secular theory with corrections arising from the 2:1 Jupiter-Saturn mean 
motion resonance. 
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Fig. 4. — Comparison between the numerical solution of the averaged equations (equations |9]- 
[TOl) and full numerical integrations of test particles at a = 2.3 AU. The dashed lines represent 
the envelope of the predicted final eccentricity, equation ( 1321) . The values of A given are in 
the canonical unit system described in ^ Each panel labeled a-d plots both an analytical 
theory result and a numerical integration result for each of the four test cases labeled a- 
d in section 13.21 The integrations were performed with Saturn starting at 8.5 AU and 
migrating outward linearly, while Jupiter remained fixed at 5.2 AU. Jupiter and Saturn had 
their current eccentricities but zero inclination. The thirty test particles in each numerical 
simulation were placed at 2.3 AU with zero inclination, but with longitudes of perihelion 
spaced 12° apart. Time zero is defined as the time when the resonance reached 2.3 AU. 
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Fig. 5. — Proper eccentricity distribution of the 931 observed asteroids with absolute mag- 
nitude H < 10.8, excluding members of coUisional families. Th e proper elements were taken 
fro m the AstDys online data service (IKnezevic fc Milanill2003l ). Family members identified 
by lNesvorny et al.l ( l2006l ) were excluded. The solid lines are the best fit Gaussian distribution 
to the observational data. The dashed line is the best fit double- Gaussian distribution. 
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Fig. 6. — The final eccentricity of asteroids as a function of asteroid semimajor axis and 
eccentricity for three different migration rates of Saturn, estimated from equation fj32|) . As- 
teroids swept by the resonance can have a range of final eccentricities depending on their 
apsidal phase, Wi. The outermost shaded region demarcates the range of final eccentricities 
for asteroids with an initial eccentricity Cj = 0.4. The innermost shaded region demarcates 
the range of final eccentricities for asteroids with an initial eccentricity = 0.2. The solid 
line at the center of the shaded regions is the final eccentricity for an asteroid with an initial 
eccentricity ej = 0. 
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Fig. 7. — The effects of the sweeping z/g resonance on an ensemble of fictitious asteroids 
with semimajor axes 2.1-2.8 AU and a uniform distribution of pericenter longitudes, (a) 
Initial distribution of eccentricities, with mean 0.05 and standard deviation 0.01 (the "cold 
belt" solution), (b) The final distribution of eccentricities after z/g sweeping, estimated with 
the analytical theory (equation ( I3T1) ). for dg = 4.0 AU My~^. The two peaks in the final 
eccentricity distribution are at approximately the same values as the observed peaks in the 
main asteroid belt eccentricity distribution shown in Fig. [51 (c) Initial distribution of ec- 
centricities, with mean 0.40 and standard deviation 0.1 (the "hot belt" solution). Asteroids 
with eccentricities greater than the Mars-crossing value were discarded, (d) The final dis- 
tribution of eccentricities after sweeping for dg = 0.8 AU My~^. The final distribution in 
d is depleted by a factor of 2.3 from the initial distribution in c. The ordinates in the four 
panels are not to the same scale. 



